Dramatic role of critical current anisotropy on flux avalanches in MgB 2 films 
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Anisotropic penetration of magnetic flux in MgB2 films grown on vicinal sapphire substrates is 
investigated using magneto-optical imaging. Regular penetration above 10 K proceeds more easily 
along the substrate surface steps, anisotropy of the critical current being 6%. At lower temperatures 
the penetration occurs via abrupt dendritic avalanches that preferentially propagate perpendicular 
to the surface steps. This inverse anisotropy in the penetration pattern becomes dramatic very 
close to 10 K where all flux avalanches propagate in the strongest-pinning direction. The observed 
behavior is fully explained using a thermomagnetic model of the dendritic instability. 



Above the lower critical field H c \ magnetic flux pene- 
trates type-II superconductors in the form of quantized 
vortices. Usually, this penetration takes place as a reg- 
ular flow leading to a smooth flux front [1]. However, 
in some superconducting films, e.g. of MgB2 a dramatic 
flux instability can occur resulting in an abrupt forma- 
tion of magnetic dendrites [2] . The flux instability is be- 
lieved to arise because the motion of flux lines dissipates 
energy leading to local heating, which results in a reduc- 
tion of pinning and an enhancement of flux line mobility. 
This constitutes a positive feedback loop which acts as 
driving force for the avalanches [3,4]. Understanding and 
controlling the avalanches is important not only from the 
fundamental point of view, but also for applications since 
they represent huge magnetic noise and reduce the effec- 
tive critical current density. 

Recently, numerous experimental investigations have 
demonstrated that the dendritic instability is sensitive 
to external conditions. The instability was found to dis- 
appear above threshold values of temperature [2] , applied 
magnetic field [5], and also for sufficiently small sample 
dimensions [6]. Also covering MgB2 films with a layer 
of gold or aluminium was seen to suppress the dendrite 
formation [7] or change their propagation direction [8] . 

To what extent the dendritic instability is sensitive 
to the internal structure and flux pinning properties of 
the superconducting film is an important and yet open 
question. A recent thermomagnetic model predicts that 
pinning strength represented by the critical current den- 
sity j c controls the onset of the dendritic instability 
[9,10]. In this paper we show that MgB2 films with 
slightly anisotropic j c can exhibit large anisotropy in 
the avalanche activity, in full agreement with the model. 
The result gives direct evidence that the microstructure 
of chemically homogeneous films grown on vicinal sub- 
strates defines a direction for dendritic flux propagation, 
which surprisingly turns out to be perpendicular to the 
lowest pinning direction. 



Substrate-induced modifications of the pinning prop- 
erties have successfully been obtained previously in high- 
temperature superconducting films grown on vicinal sub- 
strate surfaces [11,12], and on artificially patterned sub- 
strates [13,14]. To create anisotropic pinning in films 
of MgB 2 we used r-cut AI2O3 substrates slightly miscut 
relative to the [1-10 2] direction. The substrate then 
possesses a surface step structure which can be char- 
acterized by two angles AO and A(j>, according to the 
sketch in Fig. 1. Shown in the graph is the result of x-ray 
diffraction measurements obtained from surface reflection 
of a substrate mounted parallel to the 0-axis of a 4-circle 
diffractometer. The deviation 8 of the direction of the 
[1-10 2] reflection from the bisecting position is plotted 
as function of the angle <j>. From the curve it follows that 
the substrate misorientation angles are AO/2 = 0.81° and 
Acj) — 12°. With a step height of one unit cell a of AI2O3, 
the step distance becomes d = a/ tan (AO/2) m 27 nm. 
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FIG. 1. (Left) X-ray diffraction data for a sapphire sub- 
strate slightly miscut with respect to the [1-10 2] direction. 
(Right) Sketch of the square substrate with surface steps 
defining the angles A8 and A(j> which are determined from 
the plotted data. 
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FIG. 2. Magneto-optical image of the remanent state in an 
anisotropic 5x5 mm 2 MgE>2 film at 12 K. The angle a in- 
dicates the orientation of the discontinuity line which allows 
determination of the critical current density ratio j c i/jc2- 

The MgB2 films were grown by sequential deposition 
of magnesium and boron layers using a conventional elec- 
tron beam evaporation and a subsequent annealing pro- 
cess [15,16]. The films selected for this study had a lateral 
size of 5 x 5 mm 2 and a thickness of 200 nm. The mag- 
netic flux penetration into these films was investigated 
using the magneto-optical (MO) Faraday effect allowing 
direct visualization of the flux distribution [1] . In our im- 
ages the brightness represents the magnitude of the flux 
density. 

To measure the anisotropy of the critical current den- 
sity the sample was first cooled to 12 K, and then ex- 
posed to a large magnetic field which was subsequently 
removed. In this remanent state, seen in Fig. 2, the flux 
distribution develops a set of bright discontinuity lines (d- 
lines), where the current sharply bends to adapt to the 
shape of the sample [17]. In an isotropic film the angle 
between the d-lines and the edge would be 45° and the d- 
lines would meet in the center of the square. In our sam- 
ples this is not the case, as seen directly by the presence 
of a horizontal bright line segment in the central part of 
the square. Evidently, the current flowing parallel to the 
vertical edges gives stronger shielding than the current 
flowing horizontally. Assuming Bean-model behavior, 
the two critical current densities j c \ and j C 2 are related 
to the d-line angle a, see the figure, by j c i/jc2 — tana. 
Measurement of the angle gives tana = 1.06, i.e., a 6% 
anisotropy. Comparing with the direction of the steps on 
the substrate surface, see Fig. 1, the larger value of j c i 
represents a pinning force enhancement for flux moving 
perpendicular to the surface steps. 

Shown in Fig. 3 is the flux distribution after applying 
a 8 mT magnetic field to the MgB2 film initially zero- 
field-cooled to the slightly lower temperatures. Here one 
should expect the flux to penetrate in the form of den- 
dritic avalanches [2], which indeed is what we observe. 
At 8 K the avalanches are seen to enter the film evenly 



from all the edges. At 10 K, however, the penetration 
pattern suddenly becomes strongly anisotropic; flux den- 
drites develop from the vertical edges only, whereas from 
the top and bottom edges the penetration is smooth. In 
other words, the flux dendrites propagate predominantly 
in the direction perpendicular to the substrate surface 
steps, i.e., the direction of the stronger flux pinning. It 
is indeed highly counterintuitive that the flux prefers to 
move perpendicular rather than parallel to the linear de- 
fects introduced in the sample. 




FIG. 3. Flux penetration into the same MgB2 film as in 
Fig. 2, but at lower temperatures and an applied field of 
16 mT. At 8 K the flux pattern is dendritic and isotropic, 
while at 10 K dendritic avalanches are formed only along one 
direction. 



A quite similar behavior was observed for the penetra- 
tion behavior during field descent, see Fig. 4. The images 
show the remanent state after applying a maximum field 
of [1qH = 0.2 T. Again at 8 K, the penetration is isotropic 
only now the dendrites appear dark since they consist of 
antiflux surrounded by annihilation zones. At 10 K the 
dendrites propagate also now along the horizontal direc- 
tion dictated by the microstructural anisotropy. 
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T = 8 K 



T = 10 K 



FIG. 4. Remanent states of the anisotropic MgE>2 film. The 
flux pattern is isotropic at 8 K, but strongly anisotropic at 
10 K, similar to the virgin penetration behavior in Fig. 3. 

All these experiments show the remarkable fact that 
an anisotropy in j c as small as 6% can lead to a dra- 
matic anisotropy in the avalanche activity near 10 K, 
and only very near this temperature. This behavior can 
be understood using results of a recent model of the den- 
dritic instability developed in Refs. [9,10]. The model is 
based on stability analysis of the thermal diffusion and 
Maxwell equations in a long and thin superconducting 
strip thermally coupled to a substrate. Initially, the strip 
is placed in an increasing perpendicular magnetic field, 
and a Bean-like critical state is formed in the flux pen- 
etrated region. This state can become unstable with re- 
spect to perturbations in the magnetic field and temper- 
ature, and under some conditions a fastest growing per- 
turbation has a non-zero wave vector along the film edge. 
This means that an instability will develop in the form 
of narrow fingers perpendicular to the edge - a scenario 
closely resembling the observed dendritic flux behavior. 
The threshold flux penetration depth, £*, when the su- 
perconducting strip first becomes unstable, is given by 
Eq. (25) of Ref. [9], which can be written as 



Here, d is the film thickness, T* = -(din j c /dT)-\ 
k is the thermal conductivity, and ho is the coefficient 
of heat transfer from the superconducting film to the 
substrate. The parameter n characterizes the nonlinear- 
ity of the current- voltage curve of the superconductor, 
n = dlnE/dlnj > 1. 

The threshold field, H t h, when the first dendrite in- 
vades the film is obtained by combining (1) with the Bean 
model expression, 
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which connects the the applied field with the flux pene- 
tration depth £ for a long thin strip of width 2w [18]. 
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Critical current density, y c (10 10 A/m 2 ) 

FIG. 5. Threshold field for the dendritic instability as a 
function of the critical current density, calculated using (1) 
and (2). Two dashed lines show the values of j c in an MgB2 
film at 10 K for the directions along and across the substrate 
steps. Large difference in the corresponding H t h's implies a 
strong anisotropy of the dendritic avalanches specifically at 
T= 10 K. 

Fig. 5 shows the dependence of H t h on j c calculated 
using Eqs. (1) and (2) assuming kT*/E = 140 A and 
h T*/nE = 9230 A/m (which can mean e.g. T* = 10 K, 
E = 10 mV/m, k = 0.14 W/Km, n = 30, and h = 
280 W/Km 2 ). The graph shows that for high critical 
currents the threshold field is weakly dependent on j c . 
However, when j c becomes smaller, the expression in the 
brackets of Eq. (1) approaches zero, and Hth diverges. 
When j c drops below a certain value, the dendritic insta- 
bility is absent no matter how large field is applied. 

The critical current density at 10 K can be estimated 
using Eq. (2) from the flux penetration depth just 
before the avalanche behavior sets in. We then find 
jc2 = 1 x 10 11 A/m 2 , and j c i larger by 6%. These critical 
current values found from MO images agree with results 
obtained by SQUID measurements [16]. The two val- 
ues for j c are indicated by two dashed lines in Fig. 5. 
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Despite the small difference between j c \ and j C 2, the cor- 
responding threshold fields H t hi and H t hi differ signifi- 
cantly. When the increasing applied field reaches the low- 
est of the two threshold fields, H t hi, the dendrites should 
appear from the sides where j c is the highest. This is 
exactly what one can see from Fig. 3 (bottom). Interest- 
ingly, this penetration pattern dominated by dendrites 
shows an inverted anisotropy compared to the critical 
state pattern in Fig. 2. 

As the applied field is further increased and reaches 
Hth2, one could expect that dendrites appear also from 
the top and bottom edges. However, this does not hap- 
pen experimentally, as can be seen from the 10 K images. 
The reason probably is the dramatic disturbance in the 
current flow created by the dendritic structures formed 
at the smaller fields. They fill almost the whole interior 
of the film, thus lowering the edge fields substantially, 
which prevents invasion of new dendrites. 

At a higher temperature, the critical current density 
becomes smaller, i.e., we "move" to the left along the 
x axis in Fig. 5, and enter the stable region. This is 
again in full agreement with experiment, where we do 
not observe dendritic avalanches above 10 K. Conversely, 
by lowering the temperature we move to the right, where 
H t h is almost independent of j c . Hence, the dendritic 
avalanches should penetrate evenly from all four sides of 
the square film. This is exactly what we find at 8 K, see 
Fig. 4 (top), and at lower T. [19] 

Note, that the H t h{j c ) curve is slowly rising at large j c . 
It means that at very low T one should expect the oppo- 
site kind of anisotropy in the dendritic growth. Namely, 
the dendrites should first appear from the sides where 
the regular flux penetration is deeper. However, in or- 
der to observe this effect one needs a sample where the 
anisotropy in j c is much stronger than the 6% in the 
present case. 

The analysis above shows that the dendritic penetra- 
tion has a pronounced anisotropy only in the narrow in- 
terval of temperatures corresponding to the steep region 
at the H t h(jc) curve. For other superconducting samples 
this region will not necessarily be located at 10 K, but 
it will always be near the threshold temperature above 
which the instability disappears. 

In conclusion, we have investigated the magnetic flux 
penetration into MgB2 films grown on vicinal sapphire 
substrates. The resulting microstructure giving a slight 
anisotropy in j c leads to a large anisotropy in the flux 
patterns, but only in a narrow temperature interval near 
10 K. Here, the dendritic avalanches propagate only in 
the direction perpendicular to the substrate surface steps, 
i.e., the flux moves along the direction of largest pinning. 
The surprising behavior is fully explained by a thermo- 
magnetic model for the dendritic instability. 
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In principle, one should also take into account that k and 
ho are expected to increase with temperature. It can be 
easily verified that for T > 10 R this shifts the Hth(jc) 
curve to the right, while for T < 10 R - to the left. This 
produces the same relative shifts as the T-dependence 
of j c , and hence does not change the conclusions of the 
previous paragraph. 
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